A UNIFIED APPROACH TO STOCHASTIC EVOLUTION 
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Abstract. We study stochastic evolution equations driven by Gaussian noise. The 
key features of the model are that the operators in the deterministic and stochastic 
parts can have the same order and the noise can be time-only, space-only, or space- 
time. Even the simplest equations of this kind do not have a square-integrable 
solution and must be solved in special weighted spaces. We demonstrate that the 
Cameron-Martin version of the Wiener chaos decomposition leads to natural weights 
and a natural replacement of the square integrability condition. 
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1. Introduction 

Let be a standard Brownian motion, uq € L2(M), a deterministic function, and a, 
a real number. It is well known that the Ito equation 

u{t, x) = uo{x) + / Uxxis,x)ds+ / au{s,x)dW{s) (1.1) 



has a solution u{t, x) and 



/ Eu^{t,x)dx < oo (1.2) 
Jr 



for every t > and all cr G M; see, for example, [10] or simply consider the Fourier 
transform of the equation. On the other hand, the solution of equation 

u{t, x) = uo{x) + / Uxxis,x)ds+ / aux{s, x)dW{s) (1.3) 
Jo Jo 

satisfies (11. 2p for an arbitrary uq E -^^2(1^) if and only if |cr| < while equation 

u{t, x) = Uo{x) + / Uxx{s,x)ds+ / auxxi-s, x)dW{s) (1.4) 
Jo Jo 

hardly ever has a solution satisfying (11.21) . Still, because of the Ito integral, all three 
equations share a common feature, namely, that the mean evolution is described by 
the heat equation. One of the objectives of this work is to develop a unified approach 
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that would cover all three equations above without any special restrictions on the 
initial condition or the amplitude of noise. 

The situation becomes even more complicated when the noise has a spacial compo- 
nent. As an illustration, take a standard random variable ^ and consider the equation 



u{t) = 1 + / u{s)ds + / u{s)^ds. 
Jo Jo 



The solution is clearly u{t) = e*e^*. The unperturbed solution is e*, so that Ku{t) ^ e*. 

It turns out that, to preserve the mean dynamics, the usual product must be 
replaced by the Wick product uo^. In particular, the solution of 

u{t) = \^ f u{s)ds+ j u{s)o^ds (1.5) 
Jo Jo 

is u{t) = e'e^*^'-*^/^-' and satisfies E,u(t) = e*. For partial differential equations of this 
type, such as 

u{t, x) = uo{x) + / Uxxis,x)ds+ / Ux{s, x) o $^ds, (1.6) 
Jo Jo 

inequality (II. 2p . without additional assumptions on uq, holds only for sufficiently small 
t. Still, the same framework that unifies equations (11. ip . (11.31) . (II. 4p . also covers (ll.5p . 
(11.61) . and a wide class of other equations with space and space-time Gaussian noise. 

There are three main components in the general construction presented in this paper: 

(1) Use of the Skorokhod integral 6 for all types of random perturbation. 

(2) Separation of the time as the variable in the equation from the possible time 
evolution of the noise. 

(3) Replacement of the space L2(fl;X) of square- integrable X- valued processes 
with a space {C)Q^r{X) of generalized (and not necessarily adapted) random 
elements as the space for both the input data and the solution. 

Section [2] presents the necessary definitions and constructions related to the Sko- 
rokhod integral with respect to an abstract Gaussian white noise and Section [3] in- 
troduces the spaces (>C)Q,r- Section S] presents the unified treatment of stochastic 
evolution equations. While the main results. Theorem 14.71 cannot be fully under- 
stood without all the preliminaries, there is a certain analogy with the corresponding 
result for the deterministic equations. Namely, let (V, H, V) be a normal triple of 
Hilbert spaces (defined precisely in Section H]) and let A(t), t G [0,T] be a family of 
bounded linear operators from V to V. It is known from deterministic analysis, that 
if the family of operators A(t) is uniformly elliptic, then the initial value problem 
for the equation u{t) = uo + (A(s)m(s) + f{s))ds is well-posed in the normal triple 
{V,H,V')- 

Now consider the initial value problem for the stochastic equation 

u{t) = uo+ [ {A{s)u{s) + f{s))ds+ [ 6{M{s)u{s))ds. (1.7) 
Jo Jo 
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The statement of Theorem 14.71 is essentially as follows: if the family of opera- 
tors A(t) is uniformly elliptic and M(t) is a family of continuous operators from 
V to V, then the initial value problem (11. 7p is well posed in the normal triple 
{{C)Q r{V), {C)Q^r{H), {C)Q^r(y'))- The reader should not miss Example 14.41 illus- 
trating how many familiar equations are particular cases of (11. 7p . 

2. Skorokhod Integral With Respegt to Gaussian White Noise 

Let F = {Q, JF, P) be a complete probability space, and U, a real separable Hilbert 
space with inner product (■, ■)u- On F, consider a zero-mean Gaussian family 

so that 

It suffices, for our purposes, to assume that J-' is the a-algebra generated by 
W{h), hElA. Given a real separable Hilbert space X, we denote by L2(F;X) the 
Hilbert space of square- integrable JF- measurable X- valued random elements /. In 
particular, 

(/,^?)L(F;X):=E(/,^7)^. 
When X = M, we write L2(F) instead of L2(F; M). 

Definition 2.1. A formal series 

= J]iy(ufc)ufc, (2.1) 

fc>l 

where {u^, >1} is a complete orthonormal basis in lA, is called Gaussian white 
noise on U. 

Example 2.2. IfU = L2{{0,T)), thenW{h) = h{s)dw{s), where w{t) = W{x,) is 
the standard Brownian motion; Xt i^ the characteristic function of the interval [0,t]. 

Given an orthonormal basis il = {uk,k > 1} in U, define a collection {^k,k > 1} 
of independent standard Gaussian random variables by = W{uk). Let JT" be the 
collection of multi-indices a with a = {ai,a2, ■ ■ ■) so that each is a non-negative 
integer and |a| := J2k>i < oo. For a, (3 & J, we define 

a + /5 = («! + 02 + /32, • • •)) a! = ]^afc!- 

fe>i 

By (0) we denote the multi-index with all zeroes. By Si we denote the multi-index a 
with ttj = 1 and aj = for j ^ i. With this notation, nSi is the multi-index a with 
a, = n and aj = for j ^ i. 

Define the collection of random variables S = {C,a, ol G as follows: 
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where 



H^{x) = (-l)"e^V2iL,-.V2 (2.3) 



is Hermite polynomial of order n. 

Theorem 2.3 (Cameron and Martin |T]). The collection S = {C,a, ol G JT"} is an 

orthonormal basis in L2{¥;X): if t] E L2{¥;X) and r]a = ^(v^a), then rj = f]a^a 

and E\t]\^ = EaejVl 

Denote by D the Malliavin derivative on L2{¥) (see e.g. [9j). In particular, if F : 
— s> M is a smooth function and eU, i = 1, . . . N, then 

^ dF 

BFiWih), . . . WihN)) = J2 ^(^(^i)' • • • ' ^(hN))h, e L2i¥-U)- (2.4) 

i=l * 

It is known [9] that the domain D^'^(F) of the operator D is a dense linear subspace 
of L2(F). 

The adjoint of the Malliavin derivative on L2{¥) is the Skorokhod integral and is 
traditionally denoted by 5 [9J. The operator 6 extends to a subspace ©-'^■^(F; X ®U) 
of L2(F; X®U) for a Hilbert space X: given / e Di'2(F; X ®U), 6{f) is the unique 
element of L2{¥; X) with the property 

E{ip6{f))=Eif,B^)u (2.5) 

for all ip e ©^■2(F). 

We will now derive the expressions for the operators D and <5 in the basis H. To begin, 
let us compute D(^a). 

Proposition 2.4. For each a E J ^ 

D(U = 5Zv^a-.,Ufc. (2.6) 

Proof. The result follows by direct computation using the property (12. 4p of the 
Malliavin derivative and the relation H!^{x) = nHn-i{x) for the Hermite polynomials 
(cf. Chapter 1]). □ 

Remark 2.5. The set J is not closed under substraction and the expression a — 
is undefined if ak = 0. Everywhere in this paper when undefined expressions of this 
type appear, we use the following convention: if ak = 0, then ^Jakia-£k — 0- 

Proposition 2.6. For ^q, G S, /i G X, and Ufe G il, 



S{ia h®\Xk) = h Vftfc + 1 ia+sk- (2-7) 
Proof. It is enough to verify (12.51) with f = h®Ukia and ip = ^/s, where h E X. By 



Eif,-Dp)u = yMhnUp- 



^/ock + 1 h, li a = 13 - Ek, 
0, if a 7^ /? - Ek- 
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In other words, 

for all p Ej. □ 

Remark 2.7. The operator dD is linear and unbounded on L2{¥); it follows from 
Propositions \2.J\ and \2.6\ that the random variables are eig en functions of this op- 
erator: 

5(D(a)) = lal^a. (2.8) 

To give an alternative characterization of the operator 6, we define a new operation 
on the elements of S. 

Definition 2.8. For from S, define the Wick product 



e«oe.-W(^^!^)w (2.9) 

In particular, taking in (12.81) a = ksi and j3 = nSi, and using (12. 2p . we get 

Hk{^,)oK^{^,)=Hk+n{^,). (2.10) 

An immediate consequence of Proposition 12.61 and Definition 12.81 is the following 
identity: 

S{^ah ® Ufc) = h^^ o ^fc, he X. (2.11) 
More generally, we define the operation o on formal series: 



for fa e X, ga e 



zs m 



Theorem 2.9. // / = Efc>i/fc ®^k, fk = Y^a^j f^a^a e 7^L2(F;X), and f 
the domain of 5, then 

S{f) = J2fk0^k, (2.12) 

fc>i 

and 

(^(/))a = $^v^/M-e.- (2.13) 

k>l 

Proof. By linearity and (12. lip . 

^(/) = X] Si^afk,a ® Ufc) = ^ ^ fk,a^a <> ^k = ^ fk <> 6, 

fe>l aej k>l at^J k>l 

which is fl^TT^ . On the other hand, by fl^ . 

k>l ai^J k>l a&J 

and fl213|) follows. □ 
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To proceed, we need a description of a multi-index a with \a\ = n > Q using its 
characteristic set K^, that is, an ordered n-tuple = {ki, . . . , where /ci < ^2 < 
. . . < kn characterize the locations and the values of the non-zero elements of a. More 
precisely, ki is the index of the first non-zero element of a, followed by max (0, — 1) 
of entries with the same value. The next entry after that is the index of the second 
non-zero element of a, followed by max (0, — 1) of entries with the same value, 
and so on. For example, if n = 7 and a = (1, 0, 2, 0, 0, 1, 0, 3, 0, . . .), then the non-zero 
elements of a are ai = 1, = 2, ag = 1, = 3. As a result, Ka = {1, 3, 3, 6, 8, 8, 8}, 
that is, ki = 1, k2 = k^ = 3, k^ = Q^k^ = k^ = kj = 8. Note also that, for every 
sequence (6^, A; > 1) of positive numbers, 

nr=&n (2.14) 

fc>l i&Ka 

which can serve as an equivalent definition of K^. 

Using the notion of the characteristic set, we now state the following analog of the well- 
known result of Ito [3J connecting multiple Wiener integrals and Hermite polynomials. 

Proposition 2.10. Let a E J he a multi-index with \a\ = n > 1 and characteristic 
set Ka = {ki, . . . , kn}. Then 

= ^^-^^^-t:--^^^" . (2.15) 

va! 

Proof. This follows from (12.21) and (I2.10p . because by (I2.10p . for every i and k, 

k times 

□ 

By induction, we define the operator 5®", n > 1 on the space I/2(F; X Cg) W^"), where 
W®" is the symmetric tensor power of U. Then relation (12. lip leads to an alternative 
form of (|2T5|) : 

e„ = — ^5^|"|(i5„), (2.16) 

\a\' V a! 

where 

with summation in (12.171) taken over all permutations V of the set {1, . . . ,n}. This 
leads to the following generalization of the multiple Wiener integral expansion for the 
elements of L2{¥;X). 

Theorem 2.11. Ifrj G L2{¥; X), then there is a unique collection of the deterministic 
elements rjk, k > 0, with rjo = Mrj G X and rjk E X ^ U®^ , k > 1, with the properties 

^ = ^0 + E ^ '^^'(^'^)' ^ll^ll^ = ll^olll + 5^ ^11 WVkWu^. fx- (2.18) 

*;>! ■ k>l 
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Proof. Using f l2.16p and Theorem 12.31 



I — 1 \ V "-t- / 



V = Vo 

k=l "" \a\=k 

and we get the first equahty in fl2.18p with 
The second equahty in ( ]2.18p now foUows from 

k>l \a\=k 

because, by direct computation, 
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\EX^„ = k\a\. (2.19) 



□ 



3. The Spaces 

Definition 3.1. A sequence 

Q = {qk, k>l} 
is called a weight sequence if qk > 1 for all k>l. 

For a ^ J and a real number r we write 

fc>i 

Given a weight sequence Q, denote by Ag the foHowing self-adjoint operator on U: 

AgUfc = gfcUfc, k>\. 

Then, for every r G M, the operator Ag is defined. The domain of every Ag contains 
finite linear combinations of and is therefore dense in U. For / in the domain of 
Aq define the norm 

II/IIq.. = I|Aq/II^- (3.1) 

The operator Ag extends to every tensor product X ® W®"; we will keep the same 
notation for this extension and, in the case X = M, for the corresponding norm (13. ip . 
Denote by (£)'3''-(F) the dual space of (£)q^j.(1F) relative to the inner product in L2(F). 
If r/ e (/:)q,^(F;X) and ( e (/:)'5'''(F), then the duality 

U 0) 

is defined and is an element of X. 

Definition 3.2. The space (£)Q_r.(F;X) is the closure of the set of random elements 
of the form 

^ 1 - 
k=i 
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where tjq E X , r]k ^ X ®U®^, and each rjk in the domain of Aq, with respect to the 



norm 

N 



\\v\\l,r;X = WVoWx + (^11 \MQ,r \\x- (3-2) 

Theorem 3.3. A formal series rj = YlaejVaC,a is an element of {C)Q^r(^] X) if and 
only if 

E 11^-11^ <°^- (3-3) 

The left-hand side of (13. 3p . if finite, is equal to \\ri\\Q^r;X- 

Proof. By orthonormality of C,a, is enough to consider r] = fC,a, / € X. Then (12.181) 
implies 

because the definition of Ea imphes 

ll^a|||,. = g'^1«|!«!. 

□ 

Corollary 3.4. A formal series ( = Xlaej Ca^a is an element of the dual space 
(/:)Q'^(F) if and only if 

In this case, 

For h eU define the stochastic exponential 

S, = exp(w{h)-^\\h\\iy (3.4) 

Proposition 3.5. £h G (£)'^'''(F) if and only of ||/?.||q -r < 1- 

Proof. We apply Corollary 13. 4[ li h = J2k>i ^k'^i-k, hk G M, then by direct computation 
using the generating function of the Hermite polynomials, 

aej ^ fc>l 

Then, by the multinomial expansion. 



aej ' n>0 \k>l J n>0 



and the result follows. □ 
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In applications to evolutions equations, the number r in the solution space (£)(g,r;X 
is usually negative, as we want to help the series in (13.31) to converge. As a result, the 
following definition is appropriate. 

Definition 3.6. We say the the element hofUis sufficiently small relative 
to the sequence Q if there exists a positive number r such that 

\\h\\r,Q < 1- 

4. The Evolution Equation and Main Result 
Definition 4.1. The triple {V, H, V) of Hilbert spaces is called normal if and only 

(1) V "-^ H --^ V and both embeddings V ^ H and H ^ V are dense and 
continuous; 

(2) The space V is the dual of V relative to the inner product in H; 

(3) There exists a constant C > such that |(/, < C*||f ||y ||/|| y' for all v E V 
and f E H. 

For example, the Sobolev spaces {H^'^^ (R'^) , H^iM.'^) , H^f (R'^)) , 7 > 0, £ G M, form 
a normal triple. 

Denote by {v',v), v' e V , v eV, the duality between V and V relative to the inner 
product in H. The properties of the normal triple imply that \ {v',v)\ < C||f ||y , 
and, a v' E H and v eV, then {v',v) = {v',v)h- 

We will also use the following notation: 

V{T)=L,{{0,T)-V), n{T)=L2mTy,H), V'(T) = L2((0, T); \/')- (4-1) 

Given a normal triple {V,H,V'), let A{t) : V ^ V and M{t) : V V ^ U he 
bounded linear operators for every t G [0, T]. 

Definition 4.2. The solution of the stochastic evolution equation 

u{t) = u^ + (a(s)m(s) + f{s) + 5{yi{s)u{s))^ds, 0<t<T, (4.2) 

with f G Ur('^)Q>»-(^! "l^'(^)) ^^'^ ^0 £ [Jr{'^)Q,r(^', H) , is an object with the following 
properties: 

(1) There exists a weight sequence Q' and a real number r' such that 

(2) For every h E U that is sufficiently small relative to the sequence Q' , the 
function 

Uhit) = {{u{t),Sh)) 

satisfies 

Uhit) = Uhifi) + j [A{s)uh{s) + fh{s) + {m{s)uh{s),h)u)ds (4.3) 
in V{T). 
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Remark 4.3. The solution described by Definition \4-S\ belongs to the class of 
"variational solutions", which is quite typical for partial differential equations (see 
O m [HI [To] , etc.) Indeed, direct computations show that D£h = h£h, and then (14.31) 
is the result of a (formal) application of the duality {{u,Sh)) to both sides of (14.21) . 

Fix an orthonormal basis ii = {u^, k > 1} inU. Then, for every v & V, there exists 
a collection Vk E V, k > 1, such that 

M.V = ^ Vk (8> Ufc. 

k>l 

We therefore define the operators : V V hj setting M^w = Vk and write 

Mv = ^(Mfct;) ® Ufc. 

fc>i 

Then 

S{M{s)u{s)) = J2 ^k{s)u{s) {t) o ^fc, = W{uk), (4.4) 

k>l 

and equation (14.21) becomes 

u{t)=u°+ [ (A{s)u{s) + f{s) + y2^k{s)u{s)oW{uk))ds. (4.5) 

"^0 ^ k>l ^ 

If the noise W itself depends on time, this dependence is encoded in the operator M. 
As a result, (14.21) includes many popular evolution equations as particular cases. 

Example 4.4. (1) Let us see how the elementary ltd equation u{t) = 1+ j^u{s)dw{s), 
where w is a standard Brownian motion, is a particular case of (14. 2p . We take 
V = H = V' = R,U = L2((0, T)), so that W = w{t) = Ufe(t)^fe zs the Gaussian 

white noise in time. Next, put A = 0, 'WLk(t)u{t) = u(t)uk(t). Then d{M.k{t)u(t)) = 
u{t) oJ2k>i''^kit)^k = u{t)ow{t). The equivalence of u{t) = 1 + J^^ 6{'M.(s)u{s))ds and 

u{t) = 1 + Jq u{s)dw{s) now follows from the equality u{s)ow{s)ds = u{s)dw{s) ; 
see |21 Section 2.5]. 

(2) Equations with space-time white noise correspond to U = L2((0,T)) x L2{G), 
G C R"', so that the basis in U is naturally indexed by a pair of indices i,k: Uik = 
mi(t)hk{x). Thus, equation du{t,x) = Uxxdt + udW{t,x) is a particular case of (14.21) 
with yLik{t)u{t,x) = mi{t)hk{x)u{t, x) . 

(3) Equations with fractional noise are also covered by (14.21) . For example, consider 
u{t) = 1 + Jl^u{s)dW^{s), where is fractional Brownian motion and H > 1/2. 
With the appropriate interpretation of the stochastic integral, we have 

/ u{s)dw^{s)= / u{s) o {s)ds 
Jo Jo 

where 

k>l 
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and are suitable elements o/L2((0,T)) (see for example P, Chapter 5]. Then, in 
(USD, we takeU = L2((0,T)) and Mk{t)u{t) = m^{t)u{t). 

Let us fix an orthonormal basis il in W so that equation (14 ■2p becomes (14. 5p . With 
ik = ^{\ik), define ^q, a & J according to (12.21) . Given a Hilbert space X, every 
element rj of {C)Q^r{^; X) can be written as a formal series Y^^r]aia with t^q, G X 
satisfying (13. 3p . The following theorem provides a necessary and sufficient condition 
for the formal series Ylia'^a{t)ia to be a solution of (14. 5p . 

Theorem 4.5. Let u = Xlaej'^"^" be an element of {C)Qr^r'{^',V{T)) ■ The process 
u is a solution of equation ( [^.5] ) if and only if the functions Ua have the following 
properties: 

(1) every is an element of C ([0, T]; H)) . 

(2) the system of equalities 



-^0 \ k>l 



M„(t)=< + y 1^ A(s)m„(s) + + 2^v^Mfc(s)M^_,,(s) I (4.6) 

holds in V'(T) for all a E J . 

Proof. \ih = 'ZkhkUk, then, by ([33]) and Corollary EH u/, = J^a^j With no 
loss of generality we can assume that J^kyiilkY < Then, by a general result from 
functional analysis |4j, the mapping Uh : U ^ V(T) is analytic at zero. By direct 
computation, 

(1) if Uh satisfies (14. 3p . then 

1 d^'^^Uh 



Ur 



^\dh1^hf... 



h=0 



and each Ua satisfies (14. 6p . 
(2) if Ua satisfies (14. 6 p and Uh = Ylia'^ah", then Uh satisfies (14. 3p . 



□ 



This simple but very helpful result establishes the equivalence of the "physical" (14.50 
and the (stochastic) Fourier (14.60 forms of equation (14.20 . The system of equations 
(14.60 is often referred in the literature as the propagator of equation (14.50 . Note that 
the propagator is lower-triangular and can be solved by induction on \a\. 

To prove existence and uniqueness of solution of (14.20 . we make the following assump- 
tions about the operators A and M. 

(A): For every Uq & H and F G V'{T), there exists a unique function f/ G V that 
solves the deterministic equation 

U{t) = Uo + A{t)Uit)+F{t), (4.7) 

and there exists a constant C = C (A, T) so that 

||f/||v(T) < C{A,T){\\U4h + WFWv'iT)). (4.8) 
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In particular, the operator A generates a semi-group $ = ^t,s, t > s >0, and 

U{t) = ^tfiUo + [ <^t,sF{s)ds. 
Jo 

(M): For every v G V(T) and A; > 1, 







^t.,Mk(s)v(s)ds 







2 ll„.l|2 

v(r) 

V 



dt<Cl\\v\\', (4.9) 



with numbers Ck independent of v. 



Remark 4.6. There are various types of assumptions on the operator A that yield 
the statement of the assumption (A). In particular, (A) holds if the operator A is 
coercive in {V,H,V'): 

{A{t)v,v) + ^\\vfy < C\\v\\l 
for every v ^V, t & [0,T], where 7 > and C G M are both independent ofv,t. 

The following theorem is the main result of this paper. 

Theorem 4.7. In addition to (A) and (M), assume that, for some positive number 
r and a weight sequence Q, G (£)Q^_r(F; H) and f G (£)Q^_r(W; V'(T)). If 



T — 



< 00, 



then there exists a weight sequence Q° and a negative real number r° such that equation 
(14.21) has a unique solution u G (ii)Qo ,,o(F; V(T)) and 

\W\\'Q°,r°;V{T) ^ ■ (llw^llg^r;// + II / II Q,r;V'(T) 

The number C > depends only on Q,r,T and the operators A,M. 

Proof. The proof consists of two steps: first, we prove the result for deterministic 
functions m^, / and then use hnearity to extend the result to the general case. 

Step 1. Assume that the functions vP E H , f E V'(T) are deterministic so that 
u1 = , fh = f . By assumptions (14. 8 p and (14.91) the evolution equation 



Unit) =^' + (^A{s)uf,{s) + f{s) + Mk{s)uf,{s)hk^ (4.11) 

has a unique solution in V(T) as long as J2k>i ^1 small enough; since the equation 
is linear, the solution is an analytic function of h^. 

Next, 

u{t) = J2Mt)^c. (4.12) 
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and, by (14.61) . the coefficients Ua satisfy 



U( 



(o)(t)=M°+ / iA{s)u^o){s) + f{s))ds, \a\ = 0; 



u,^{t) = [ A{s)u,^{s)ds + [ Mfc(s)u(o)(s)ds, |a| = 1; (4.13) 
Jo Jo 

Ua{s) = / A(s)uQ,(s)(is + a/o^ / Mk{s)ua-ek{s)ds, |a| > 1. 
Jo ^ Jo 

Denote by $ = ^s,t, t > s > the semigroup generated by the operator A(t). It 
follows by induction on |q;| that 

M(o)(t) = <^t,ou^{x) + [ <!>t,sf{s)ds, \a\ = 0; 

u^^(t) = / ^t^sMk{s)u(o){s)ds, |a| = 1; 
Jo 

Mt) = ^= / •••/ ^t,sM^n^Sr.,Sr..^■■■Mi^^_^^s2,s-M^lU{0)dSl...dSn, 

V al Jo Jo Jo 

(4.14) 



\a\ = n> l,Ka = {ii, . . . ,in} 
By assumptions (14.81) and (14.91) . 



MnliT) < ^ {\\A\h + ll/llv'(T)) n^'f (4.15) 



k>l 



It is known (see, for example, [21 pages 31, 35]) that 

|a|! < a!(2N)2", ^(2N)~™ < oo if and only if r > 1; (4.16) 

where 

k>l 

Define the sequence Q° = {ql, > 1} by 

ql = 2kil + Ck). (4.17) 

Then (0 and fHlSD imply 

lkllQo,.°;V(T) < C ■ {\\uYh + ll/Wllv'(T)) for every r° < -2. 

Step 2. As in Step 1, existence and uniqueness of solution follows from unique solv- 
ability of the parabolic equation (14. lip , and it remains to establish (14.101) . 

Given v e H, F e V'(T), and 7 G JT, denote by u{t; v, F, 7) the solution of (C2|) with 
M° = / = ^^7- If «° = Eaej / = T^aej by linearity, 

M(t) = ^n(t;<,^,7). (4.18) 
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It follows from f l4.13p that Ua{t; v, F, 7) = if |a| < I7I and 



(4.19) 



^/{a + 7)! a/oI 
Using the results of Step 1, 

lk(t;<, A,7)llQ°,r-°;V(T) < (11^711/^ + IIAWIIvm) • (4-20) 

With Q° defined in ( KT7\i and r° < -2, fli?TO|) now follows from (14181) by the triangle 
and the Cauchy-Schwartz inequalities. □ 

Corollary 4.8 (A generalization of the Krylov-Veretennikov formula from [6]). As- 
sume that M° and f are deterministic and take from (14.171) . Define the sequence 
Un = Unit) , n > by induction 

Uoit)=u^o)it), f/„+i(t)= [ <^t,sS{M''{s)Ur,{s))ds, n>0, (4.21) 

where M° = (g^Mi, q°M2, ...). Then 

Y,{qTua{t)i^ = Un{t). 

\a\=n 

Proof. By Km . 

{qTMt)= f A{qTua{s)ds + J2 f qkV^kM,{qT~''Ua-e,{s)ds. 
Jo Jo 

Therefore, 

(f/n(t))a = 5^V^ / $t-.M°(f/„_i(s))„_,,rfs. 

k>l 

By (I2I3]), 

Un{t) = [ <!>t,sS{M°{s)Un-l{s))ds 

Jo 

and the result follows. □ 

Remark 4.9. The main goal of Theorem \4.7\ is universality: to cover the largest 
possible class of equations. It is therefore inevitable, and natural, that, for many 
particular equations, there are much better regularity results than (I4.10p . 
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